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Abstract 

We have designed a new experimental setup able to investigate fracture of 
soft materials at small scales. At high crack velocity, where energy is mostly 
dissipated through viscoelastic processes, we observe an increasingly large high 
strain domain in the crack tip vicinity. Taking advantage of our ability to 
determine where linear elasticity breaks down, we derive a simple prediction for 
the evolution of the energy release rate with the crack velocity. 
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1. Introduction 

Soft materials are of great interest both for industrial and fundamental rea¬ 
sons. Materials with shear moduli in the range 0.1-100 kPa can be found indeed 
in very different applications, from tissue engineering to food industry or geo¬ 
physics. While ubiquitous, their mechanical properties are poorly understood, 
especially when related to fracture. First, it is an experimental challenge to 
perform controlled fracture experiments on such materials, which are difficult 
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to grip [? ] and which, sometimes, flow at a macroscopic scale [? ]. Second, 
the fracture mechanics of these materials is made particularly complex by the 
interplay of large deformations [? ] and viscoelastic processes [? ]. 

Fracture of soft materials has been widely studied both theoretically [? ? 
? ] and experimentally [?????]. Large efforts were made to measure 
the rate dependency of the energy release rate [? ? ? ] and to model it using 
either viscoelastic models [? ? ] or microscopic damage models [? ]. However, 
in very deformable materials, microscopic cohesive zone models deal with what 
happens in the very non linear vicinity of the crack tip [? ? ] • Because non linear 
processes will change the way the elastic energy is conveyed toward the crack tip 
and then dissipated in its vicinity, it is crucial to characterize crack propagation 
at both macroscopic and microscopic scales. Biopolymer physical gels are good 
candidates to gain a deeper understanding of the fracture properties of soft 
materials. Indeed, linear and nonlinear mechanics of these gels have been widely 
characterized [? ? ]• 

In this Letter, we introduce a new experimental setup that was designed to 
probe the mechanical response of soft materials in the vicinity of the crack tip. 
As described in Section Mode I cracks are grown at controlled velocity V in 
a 2D sample of agar gel. In Section we analyse crack profiles when varying 
V. We show that Linear Elastic Fracture Mechanics (LEFM) breaks down at a 
distance from the crack tip which increases with V and which is proportional to 
the viscoelastic dissipated energy. Finally, Section is devoted to a discussion 
of the experimental results. 

2. Material and Methods 

Material and sample preparation. Agar gels are prepared by dissolving agar 
powder (Sigma Aldrich, average molecular weight 12{)kg.mol~^ ^ chain contour 
length A = 200nm) in hot water. A 1.5%wt solution of agar is stirred for Ih at 
95 C without any water loss. The obtained sol is then kept at 70 C and used 
within a day. The gelation mechanism is simple: at temperatures lower than 
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37 C, dissolved chains start interacting via H bonds and locally form junctions, 
leading to the formation of a weak physical network. 

Linear dynamical rheology measurements are performed at 22 C (Anton 
Paar Rheometer, 1% strain in the linear viscoelastic range, roughened cone 
and plate geometry), showing that both G' and G” are nearly constant within 
the frequency range (0.1 Hz-100 Hz), with G” being 40 times smaller than G' 
(G'=20kPa). The material being incompressible, its Young modulus is E = 
3G'=60kPa. Knowing E gives us an order of magnitude for the typical distance 
between the junctions ^ = (kT/EY^^ lOnm. The non linear elasticity of agar 
gels was investigated by Pavan & al [? ]. This work shows that agar gels are 
not neo-Hookean but would rather exhibit exponentially stiffening. 

Just before casting the gel, 1ml of sol is dyed with 2.5/il of a dispersion of 
carbon black in indian ink (inset of Fig. |^. 



Figure 1: On a chip mechanical test cell. Hele-Shaw chamber containing the sample of 
gel (in red) surrounded by fluorinated oil (in blue). Oil is sucked out of the chamber from 
points A, which results in a displacement of the gel boundaries in the y-direction. In this 
plane stress configuration, a mode I crack grows in the x-direction. 


Fracture experiments. Our new device is composed of a rigid chamber made of 
two glass plates of thickness 1mm, separated by a distance of 350/im (Fig. [^. 
It is microfabricated in a clean room, using microfluidic stickers technology [? 
]. The edges of the chamber are sealed with photocurable glue. 

It is first filled with FC3280 fluorinated oil (in blue in Fig. [^. The hot 
sol (in red in Fig. is injected in the middle of the chip through a hole B in 
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the upper plate, and cooled down to 22C for 1 hour until complete gelation. 
Before starting the fracture experiment, the gel sample dimensions are 0.35mm 
(thickness) in the z-direction, 25mm in the x-direction and 14mm in the y one 
(Fig. 0. To assure perfect slippage and prevent adhesion of the gel on all the 
walls of the chamber, the latter were coated with an acrylamide brush using a 
two-step surface reaction [? ]. Two synchronized syringe pumps connected to 
points A (Fig. suck the oil out of the chamber at a prescribed flow rate Q. 

The oil being immiscible with the gel, the interfaces between the fluid and 
the material remain sharp. Acting as a tensile machine of low stiffness, the 
oil outflow imposes a velocity of displacement proportional to Q to these two 
boundaries. However, due to the low rigidity of the tensile system, the dis¬ 
placement along the oil/gel interface is non uniform as soon as a crack starts 
propagating. 

Under deformation, a single mode I crack is nucleated at point B from a 
built-in notch of length 5mm and radius of curvature at the tip lOjam. The 
crack tip is observed under a microscope with a 4x magnification (observation 
field 9 mm^, pixel size 1.5/im). 

B is connected to a water reservoir so that the propagating crack is filled with 
water instead of air. As a consequence, our gel being mostly made of water, 
we suppress surface tension forces which might generate a spurious blunting 
of the crack tip, smoothing it at scales smaller than the elastocapillary length 
j/E = 1.2/im. 

Experimental observations. At fixed Q, after a transient acceleration, the crack 
reaches a propagation regime where the crack velocity V varies very slowly with 
time, until it starts feeling the edges of the specimen. Within the steady-state 
regime, sequences of images are captured and processed with a home-made 
image analysis routine detecting the crack lips. At this concentration, changing 
the control parameter Q from 0.1 to 1000/iL/min enables one to tune V between 
l/im.s“^ and 1cm.s“^: V being always much smaller than the Rayleigh wave 
speed Vr lm.s“^, crack propagation can be considered as quasi-static. Note 
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that for velocities lower than crack propagation becomes intermittent: 

under stress, junctions have time to unzip and the network to relax. For V higher 
than 1 cm.s~^, it becomes difficult to reach a steady state. 

For each experiment, the crack opening displacement (COD) is extracted 
from the observed shape of the crack, while the full displacement field around 
the crack tip is determined by using the CorreliQ4 Digital Image Correlation 
(Die) code [? ]. Full DIC results will be reported elsewhere [? ], where the 
dominance of mode I over mode II (Kj=50 Ku) is shown. In this plane stress 
configuration, we can thus relate directly the energy release rate S= KjjE [? ] 
with the stress intensity factor (SIF) Kj. Kj can indeed be extracted from the 
measurement of the displacement along the crack direction, behind the crack 
tip. 


3. Experimental results 
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Figure 2: Slow crack Crack opening (displacement along the crack lips) as 

a function of the distance d to the crack tip. Red line with circles: experimental results. Black 
dashed line: fit with the LEFM expression^ Inset: Picture of a low-velocity crack (white) 
in a dyed 1.5%wt agar gel. The interior of the crack is filled with water. The contour of the 
crack is superimposed in red. 
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Fig .[^shows a typical snapshot of a crack in the low velocity range {V =3/ini.s 
and its corresponding COD. A video sequence can be accessed from the link pro¬ 
vided in [? ]. This profile fits in well the LEFM theory and Williams’s prediction 
(1957) for the crack opening displacement u\ 




d 


l + ^+ ^ 


di 


d 


d2 


( 1 ) 


where d is the distance to the crack tip, and E is the gel’s Young modulus. 
Fig-i shows a fit of the experimental profile based on the first three terms 
of Williams expansion. For this low-velocity crack, di and d 2 are respectively 
900/im and 4mm, showing that, even if they are essential to capture the crack 
profile, these two extra terms are only relevant far from the crack tip. This 
shows the existence of a zone of i^/-dominance at the crack tip, i.e. a region of 
space in which the Irwin term dominates [? ]. Hence Kj can be deduced 
directly from the COD measurement. 
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Figure 3: Rapid crack (y=lcm.s“^). Crack opening (displacement along the crack lips) 
as a function of the distance d to the crack tip. Red line with circles: experimental results. 
Black dashed-dotted line: fit of the data with the LEFM prediction ^ which clearly does 
not work close to the crack tip; Blue dashed line: fit of the data with the modified Williams 
expansion!^ Inset: Picture of a high velocity crack in a dyed 1.5%wt agar gel. The interior 
of the crack is filled with water. The contour of the crack is superimposed in red. 
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However, when increasing V, the crack shape departs from the LEFM pre¬ 
diction. Fig.j^shows the COD for a high velocity crack (F=lcm.s“^). A video 
sequence showing the crack propagation can be accessed from the link given 
in [? ]. The crack profile cannot be captured by Williams expansion. Instead of 
an “infinite” slope (at the experimental resolution) at the LEFM crack tip, fast 
cracks propagating in agar gels exhibit a finite slope at the apex with u(d) oc d. 
A central point of this paper is to determine the very nature of the region close 
to the crack tip where this is observed. Is it a Dugdale Barenblatt cohesive 
zone [? ]? Is this new crack shape due to the material viscoelasticity — the 
so-called de Gennes trumpet [? ]—? Or is it an effect of large deformations 
at the crack tip [? ]? Based on our phenomenological observation, we write a 
modified Williams expansion series: 


u{d) = d 


A 

l + (|r)V2 



( 2 ) 


where d^ is a crossover length between the linear apex and the classical LEFM 
expansion. A is a measure of the finite slope at the apex. While the crack 
opening profile may look like a perfectly linear elastic one at large scale, it is 
obviously not the case at distances smaller than d*. Eq.j^fits well the COD over 
the whole velocity range. di and ^2 hardly change with crack velocity. They 
are both of order 10^/im with di < ^2 : the corresponding terms dominate at 
distances larger than 1mm. Each crack propagating at velocity V will be fully 
characterized by a couple (A(F), d*(F)). 
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Figure 4: Parameters of the modified Williams expansion as a function of crack 
velocity V. Left. Evolution of with V. Right Evolution of A with V. The blue line is a 
guide for the eye showing that A is nearly E-independent. 

Values of A and d* are shown in Fig. Except for very low velocities, for 
which it is difficult to define a linear apex, A appears to be a constant of order 
unity over the whole velocity range. Meanwhile, a continuous increase of 
with V can be observed (Fig.[^ left), appears to be a dynamical lengthscale 
that vanishes at small velocities and may reach up to 100 fim at high velocity. 
Similar results have already been observed by Bouchbinder & al [? ] in the 
context of the dynamic fracture of chemical gels. It is the first time, however, 
that such an observation is made for quasi-static crack propagation in a soft 
material. 

At distances from the crack tip which are large compared to d^, but for 
which the 2^^ term of Eq. is still negligible, the dominant contribution to 
the opening is ^ A{d'^dy^‘^. By analogy with the Irwin term in Eq. we can 
compute the SIE Ki{V): 

Ki{V) = A{V)d*{vA^E^ (3) 

The energy release rate S is computed using Eq. i and S= (see Eig.[^. 

The evolution of S as a function of V is shown in Fig. As expected for such 
soft materials where viscous processes occur, S increases with V due to viscous 
dissipation within the bulk of the material, or due to dissipative processes at the 

8 





crack tip. At low velocity V, S seems to reach a finite limit Sc below O.lJ/m^. 
It can be interpreted as the critical energy release rate at the onset of crack 
propagation. At high velocity, S ^ Sc- Increasing 1^ by 4 decades induces more 
than one order of magnitude increase of S • We note that the results obtained 
for S (V) by analyzing crack profiles in agar gels are qualitatively similar to 
those obtained by integrating force-elongation curves measured during fracture 
experiments on other biopolymer or physical gels [? ? ]• 

4. Discussion. 

A natural assumption would be that d* is the size of a process zone where 
viscoelastic dissipative processes take place. Dugdale—Barenblatt models [? ] 
consider a process zone of size £j:)b where cohesive stresses prevent the crack 
from opening. But, whatever the cohesive stress profile, asymptotic crack open¬ 
ing is predicted to show a cusp in [? ] which does not match our experi¬ 
mental observations at this scale. Besides, as G' > G” over the whole frequency 
range, the material is elastic at all time scales and there is no reason for the 
crack tip shape to result from bulk viscoelasticity. 
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Figure 5: Evolution of S with V. Black circles are the experimental results computed using 
Eq. The red dashed line is a fit corresponding to S= Sc(l + Inset. Linear 

variation of S with the dynamic lengthscale d*. 


Finally, let us assume that marks a domain of large deformations in the 
crack tip vicinity, which leads to a strong deviation from LEFM. 

From Eq.[^ at a distance r > ahead of the crack tip, i.e. within the linear 
elastic zone, the deformation scales as: 


e(r) = dru{r) 


Ki{V) 

E\^27rr 


A{V) 


( 4 ) 


A being independent of E, with A ^ 1.2, we find that e{d'^) is about 30 %, 
which is a deformation large enough to justify that LEEM does not apply in 
this zone. Indeed, Bouchbinder & al [? ] have shown that the effect of large 
deformations becomes measurable within the 10-20% range. At distances d < 
d*, the crack tip shape depends strongly on the nonlinear constitutive law of the 
material. While, for neo-Hookean materials, Geubelle & Knauss [? ] predict a 
” linear-like” parabolic crack tip which was observed in [? ], Long & al compute 
the non linear crack tip shape in an exponentially stiffening material [? ]: they 
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predict a wedge-like shape similar to what is observed in our case. 
We now discuss the ^/-dependence of S which can be written as: 


S = Sc + 9vis{V) 


( 5 ) 


where Sc is the critical energy release rate (fracture energy at vanishing crack 
velocity) and Svis(^) is the rate-dependent energy (per unit crack surface cre¬ 
ated) dissipated during steady-state propagation at velocity V. 

Following [? ], we estimate the bulk dissipation for a stationnary crack (for 
which time derivatives can be replaced by space derivatives, i.e. Vdt = dr) 
propagating at vanishing velocity V in the following way: 


9vis = ^ J {drufrdr ( 6 ) 

pr=oo 

= r]V / {drru)‘^rdr (7) 

Jr=dc 

where we have not taken into account the angular variation of u since, once 
integrated in Eqj^ it will only modify the prefactor. In [? ], 77 is an effective 
material viscosity, dc is the microscopic length scale where the linear theory 
fails, which is actually the lengthscale d'^ measured in our experiment. 

In this perturbative approach, one can use the elastic asymptotic crack dis¬ 
placement field u{r) oc {Kic/E)y/T^ and oc —Kicj Hence: 


p Sc^ _ 

9vis oc 


(8) 


where Sc = 

Computing 9vis{V) at higher velocity is a more delicate task. Indeed, at 
finite H, local dissipation at a distance r from the crack tip depends on the 
local deformation rate at pulsation uj = V/r [? ]. Moreover, the perturbative 
approach, valid when S^is Sc does not hold anymore when S ^ Sc- If we 
recall that e* = e(r = d'^{V)) is velocity-independent, then Eq. directly leads 
to: 


S(I") = 


KdiV) 

E 


2nE{e*fd^V) 


( 9 ) 
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i.e. S{V)/E oc (E. The inset of Fig. shows S versus (E. A linear relation 
9{V) = ad^ fits indeed the data with a = 2.7 lO^Pa . For large y, as S ^ 9vis, 
we have: 


d-{V) 


9c d*iV = 0 ) 

Matched asymptotics enables to eliminate d*{V) in Eq. [fusing Eq. 10, and 
we get, for high velocities: 


riV 


''\lEd*{V = 0) 

Over the whole velocity range, we thus have: 


S ~ Sc(l + 


7]V 


Ed*{V = 0) 


) = 5cil + xl^) 


( 11 ) 


( 12 ) 


where 


V = 


Ed*(V = 0) 


(13) 


Agreement between the measured fracture energy 9(V) and Eq. 12 is good 
within the explored velocity range (Fig. [^. 

The resulting value for the critical energy release rate 9c = 0.01 J/m^ is 
reasonable. Let us recall that 9c is the surface energy needed to break the 
bonds crossing the fracture plane. Given the weakness of H bonds compared 
to the covalent bonds, fracture of these physical gels is reported to occur via 
unzipping of the junctions of energy Uh [? ]• The model of Lake & Thomas 
(1967), stating that all the energy stored in a chain containing n residues is 
dissipated at rupture, 9c should scale like UhI^cl where ^ is the distance between 
the junctions and a is the size of a residue. With Uh = O.leE, ^ = lOnm and 
a = 0.3nm, we find 9c = 5.10“^ J/m^ which is in good agreement with the value 
measured thanks to our model. 

The fit also gives a value for V = 0.1 which is below the minimum 

crack velocity that we could reach. For a crack of velocity E, S is expected 
to be twice Sc- What we observe experimentally is in fact an intermittent 
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crack propagation regime. At these velocities, we suspect that a moderate 
stress relaxation of the network is sufficient to dissipate elastic energy and to 
decrease S below its critical value Sc causing the crack to grow intermittently. 

we get T] = Ed*{V = 0)/!^. When no large deformation 


From Eq. 
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is expected at the crack tip and the limitation to LEFM is thus the extent 
of the Dugdale—Barenblatt cohesive zone Idb = ^cl(^y = Sc/E 0.1/im. 
Consequently, r] is found to be EIob/V ^ 6.10^Pa.s. 

We can wonder how a material with rate independent G' and G” exhibit rate- 
dependent fracture properties. The answer lies in the rate-dependency of the 
highly non linear processes occuring in the vicinity of the crack tip. One of them 
is viscous chain pull-out, proposed by Baumberger in the case of gelatin [? ]. It 
results in an effective viscosity scaling like r]sE/^ = 20r]s where r]s = 10“^Pa.5, 
is the solvent viscosity because we consider here chain/water friction. Eor our 
agar gel, this effective viscosity is orders of magnitude too low. 

Another option is to consider that dissipation stems from the dynamics of 
agar chain segments confined in the junction zones. Indeed, reptation processes 
result in an effective viscosity r]ef f = GGrep [? ] where Trep is the chain reptation 
time. In our case, the reptation time is the time at which our biopolymer 
network relaxes. At our working concentration, [? ] found that agar networks 
relax on timescales of 10^ — 10^5. This would result in r]eff = 10^ — 10^Pa.s 
which is now far larger than our measurement and confirms that our material 
does not flow on experimental timescales. But if we now consider the material 
in the vicinity of the crack tip, stress-aided reptation of the strongly deformed 
chains may significantly decrease Trep by I or 2 orders of magnitude [? ], which 
leads to rjeff in the measured range. 

Finally, the observed increase of (P with V suggests an increase of the defor¬ 
mation to failure for increasing crack velocity and confirms the rate-dependency 
of the failure properties of agar networks (see for example Fig. 5, 6 and 7 in [? 
]). Below a velocity of order V (which is out of our experimental range), stress- 
aided partial unzipping of the junctions leads to an increased strain to failure 
cr of the network. At moderate F, cr reaches a mininum lower than 10%: 
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Fracture occurs in the linear elastic regime, and d* is thus very small. Finally, 
at higher rates, eR increases with V^ and reaches values larger than 10%. This 
is the reason why we observe a more extended zone of high strains in the crack 
tip vicinity when V is much larger than V. A quantitative understanding of 
the rate depency of polymer gel failure properties at high deformation rates still 
lacking. 

5. Conclusion. 

In our experiments, we were able to recover the energy release rate measured 
by other authors [? ? ? ] from macroscopic experiments, by analyzing the 
morphology of the crack tip. 

Because fracturing our physical gels at high velocity V causes a great deal 
of dissipation, the stress intensity factor Kj must increase with V accordingly. 
This induces large deformations in the crack tip vicinity, which translates into 
an earlier departure from LEFM (d* increases with V). A full understanding of 
the very shape of the crack at small scales is still in progess, based on the work 
of Livne & al [? ] and Long & al [? ]. 

We will take advantage of our new experimental setup to study less concen¬ 
trated gels with lower moduli, in order to observe the emergence of the cohesive 
zone. More fundamentally, the growth and possible divergence of the size of 
the process zone in the vicinity of the sol gel transition will be studied. In 
the future, these experiments can be extended to other types of liquid to solid 
transitions such as the colloidal glass transition. 
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